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Abstract 

A generalized diffusion Monte Carlo method for solving the many-body 

Schrodinger equation on curved manifolds is introduced and used to perform a 

'fixed-phase' simulation of the fractional quantum Hall effect on the Haldane 

sphere. This new method is used to study the effect of Landau level mixing 

on the u = 1/3 energy gap and the relative stability of spin-polarized and 

spin-reversed quasielectron excitations. 
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Much of the modem understanding of the fractional quantum Hall effect (FQHE) is 
based on the observation that in two dimensions the quantum statistics of identical particles 
can be changed by performing a singular gauge transformation [Q. Such a transformation 
can be used, for example, to map the equations describing an ideal two dimensional electron 
gas in a transverse magnetic field at Landau level filling factor v = where q is an odd in- 
teger, into those describing a system of 'composite bosons' moving in zero effective magnetic 
field, interacting via both Coulomb and 'Chern-Simons' interactions. This transformation, 
from fermions to bosons, is the basis of the successful Chern-Simons-Landau-Ginzburg phe- 
nomenology of the FQHE p|. Its existence also suggests the possibility of numerically 
simulating the FQHE at v = 1/q using the composite boson description. 

Recently, Ortiz, Ceperley, and Martin (OCM) [Q] introduced the 'fixed-phase' diffu- 
sion Monte Carlo (DMC) method for simulating non-time-reversal symmetric systems with 
complex-valued eigenfunctions. OCM applied this method to the FQHE ground state at 
V = 1/q, using the torus geometry, and fixing the phase of the wave function with Laugh- 
lin's trial wave function. The resulting effective bosonic problem corresponded precisely 
to the composite boson description, with the additional approximation that those terms in 
the transformed Hamiltonian leading to fluctuations of the phase of the wave function were 
ignored. This effective bosonic problem was then solved by standard DMC techniques [Q, 
and the results used to study the effect of Landau level mixing (LLM) on the FQHE ground 
state. However, OCM did not consider either excited states or geometries other than the 
torus. 

In this Letter we present the results of a 'fixed-phase' DMC study of the FQHE at 
u = 1/3 using the spherical geometry introduced by Haldane 0. In order to go from the 
torus to the sphere we introduce a generalized DMC method for solving the many-body 
Schrodinger equation on curved manifolds. One motivation for this work is that the 'Hal- 
dane sphere' is arguably the most convenient geometry for numerical study of the FQHE, 
and we believe that our generalization of the fixed-phase DMC method to this geometry 
will be useful for many future calculations. As an example of the application of this new 
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method we have calculated the effect of LLM on the FQHE transport gap, i.e., the energy 
gap for creating a fractionally charged quasielectron - quasihole pair with infinite separa- 
tion. Results have been obtained for both a spin-polarized and spin-reversed quasielectron. 
Previous calculations of the crossover magnetic field below which the transport gap is set by 
the spin- reversed excitation have ignored LLM . The present work includes these effects 
for the first time. 

An ideal two dimensional electron gas, with effective mass m*, carrier density n and 
dielectric constant e, placed in a transverse magnetic field B is characterized by three length 
scales — the effective Bohr radius as = efi^ /m*e^, the magnetic length Zq = \J'hc/eB, 
and the mean interparticle spacing d = 1/ ^Jwn. These length scales can be combined to 
form two independent dimensionless ratios, the filling factor v = and the electron 

gas parameter = d/as- The degree of LLM is characterized by the ratio of the typical 
Coulomb energy to the cyclotron energy {e^ / ed) /huOc = TsZ//2 where Uc = eB/m*c. Thus, for 
fixed Tg provides a useful measure of the importance of LLM. Because oc m* /\/B LLM 
can be increased either by decreasing B or increasing m* . For example, in two dimensional 
GaAs/AlGaAs systems with typical carrier densities, for n-type systems m* ^ 0.07 and 
~ 2 while for p-type systems m* ~ 0.38 and ~ 10. 

In the spherical geometry electrons are confined to the surface of a sphere of radius R 
with a magnetic monopole at its center. Let N denote the number of electrons and 25* 
denote the number of flux quanta piercing the surface of the sphere. The field strength is 
then B = STic/eR^ and, for v = 1/q, 2S = q{N — 1). If electron positions are given in 
stereographic coordinates, r = {x,y) = (cos 0, sin0) tan6'/2 where 6 and are the usual 
spherical angles, then the Hamiltonian is, in appropriately scaled atomic units {h = m* = 
e' = 1), 

i/ = iEfWHV, + A(r,)f + li:^=;. (1) 

where -D(rj) = (1 + rfY /AR^, rij is the chord distance on the sphere, and /? is a parameter 
which models the finite thickness of the 2DEG We work in the Wu-Yang 0] gauge for 
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which Ain) = 2R^B{-yi,Xi)/{l + rf). 

In what follows we have used three trial wave functions to implement the fixed-phase 
approximation {z = x + iy is the complex stereographic coordinate): 

(i) Ground state. — In the Wu-Yang gauge the spherical analog of the Laughlin wave 
function for z/ = 1/g is H 



^Gs=n(i+i^'cp)-^n( 



(2) 



i<j 



(a) Spin-polarized excited state. — For this state we have used the following wave function, 
constructed using Jain's composite fermion approach |T^, describing an excitation with a 
charge e/q quasielectron at the top of the sphere and a charge —e/q quasihole at the bottom 
of the sphere, 



k i<j 



1 ^ V- 1 

i Zi ... z^i^i — 



i Zat ... Zj^ l^ij^N 



(3) 



(Hi) Spin-reversed excited state. — This state is similar to t/'sp except the quasielectron 
has a reversed spin. If zi denotes the coordinate of the down spin electron then this wave 
function is 



Zi Zj J 



(4) 



k Ij^l i<j 

The fixed-phase approximation is carried out as follows. First, the relevant trial func- 
tion is written as iPt(J^) = \4't(J^)\ exp[z0T(7^)], where 71 = (ri, r2, ■ ■ ■ , tat). The 
overall phase 0t('7^) is then used to perform a singular gauge transformation, H = 
exp[— 'i0T('^)] H exp[i(j)T(Jt)] = H^ + iHi, where 



Hk = -W D(r.)(V? - A^(r.)) + \ ^ 



Hi = -\T. D{vi) (v. ■ A(r,) + A(r,) ■ V, 



(5) 
(6) 



and A(r.j) = A(ri) + Vi0T('^)- As shown by OCM, the bosonic ground state of is the 
lowest energy state with the same phase as the trial function [Q. The DMC method can 



then be applied to the imaginary time Schrodinger equation Hji'ijj{Tl,t) = — ^-0(7^, t), the 
solution of which, in the limit t oo, converges to the ground state of i^R. 

The dependence of -D(r) on position in (|I]) is a consequence of the finite curvature 
of the surface of the sphere, for which the metric tensor, in stereographic coordinates, is 
Qapiy) = D{Y)^^5ai3- Bclow we introduce our generalized DMC method for simulating the 
many-body Schrodinger equation on such a curved manifold. Note that in two dimensions 
it is always possible to choose coordinates for which gap = f{r)6aj3, i-e., to work in the 
so-called 'conformal gauge', and so the generalized DMC method introduced below can, in 
principle, be applied to any curved two-dimensional manifold, not just the sphere. 

The central modification of the DMC method required to treat the present problem is 
to replace the usual importance-sampled distribution function P(7l,t) = ■ip(Jl,t)\ip'T(JZ)\ Q 
with 

^ 1 

p{n,t) = 4j{n,t)\MT^)\ U 7^ • (7) 

i=i '^{'^i) 

This has two important consequences. First, because the differential area element on the 
sphere is dA = d'^r/D{r) the expectation value of the ground state energy is simply 

^ /p(7^,t^oo)Ez.(7^)rf7^ 

jp{n,t^ oo)dn ■ ^ ' 

Second, the differential equation satisfied by -P(7^, t) is, 

N r 



-lv,^(z^(r,)P(7^,^)) + V, ■ (D(r,)F,(7^)P(7^,^)) 



+ {EL{n) - ET)P{n,t) , (9) 



where Fj(7^) = Vjinl'j/'rl, EiiTZ) = H^\ipT^\/\tp'j^\ and is a constant which must be 
adjusted in the course of the simulation to be equal to the ground state energy. It is worth 
noting that, except for the position dependence of -D(r), (|) has the same form as the usual 
generalized diffusion equation appearing in DMC simulations 0]. 

Equation can be solved numerically by stochastically iterating the integral equation 

p(7^', t + T) = J G{n 7^', r)P(7^, t)dn , (lo) 
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using the short-time propagator ((9(r^)) 

f[EL{n) + EL{n') 



exp 



N 



^G°(7^-7^^r) 



(11) 



1=1 



where 



G°(7^^7^^r) 



2TTD{ri)T 



exp 



(r^-r,-D(r0rF,(7^))- 



(12) 



represents a diffusion and drift process. In (|T2]) both -D(ri) and FiiJZ) are evaluated at 



the 'prepoint' in the integral equation. This makes it possible to simulate (|iy) in terms of 
branching random walks by a straightforward application of the rules given in and is a 
direct consequence of having the spatial derivatives in sit all the way to the left in each 



term ||TT|. This in turn follows from our modified definition of P(lZ,t). Had we used the 
usual definition, P(7l,t), we would have obtained 'quantum corrections' in the propagator 
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The Tg dependence of the ground state energy, and the energies of the spin-polarized 
and spin-reversed excited states, have been calculated by fixing the phase with the trial 
functions ipcs, V'sp and ipSR, respectively, and solving the resulting bosonic problems using 
the generalized DMC method outlined above. Figure 1 shows the results for the ground 
state energy as a function of for (3 = and u = 1/3, compared with the variational 



Monte Carlo results of Price et al. |jT3|,[I^. This comparison provides an important test of 
our generalized DMC method — the wave functions used in the variational calculations have 
the same phase as ipQs and so must have higher energies than our fixed-phase DMC results, 
as is in fact the case. 

The spin-polarized and spin-reversed energy gaps, Asp and Asr, obtained by subtracting 
the ground state energy from the relevant excited state energies, are plotted vs. in Fig. 2. 
Results are for = 20 and are given for (3 = and the more experimentally relevant case 
P = 1.5 Iq. To reduce finite size effects we have subtracted Vq = —(6/^)^/26/2, the Coulomb 
energy of two point charges with charge ±e/g at the top and bottom of the sphere, from 
our results for the gaps [|l^]. The additional Zeeman energy of the spin-reversed excitation 



is not included in our definition of Asr, and so tlie crossover magnetic field, Be, below 
which the spin-reversed excitation has lower energy than the spin-polarized excitation, is 
Be = (Asp — AsR)/gfiB, where /ib is the Bohr magneton and g is the effective (7-factor. 
For GaAs (e ~ 13,(7 ~ 0.5) we find, for Vs = 2, Be — 14 T, while for = 10, Be — 7 
T. This reduction of Be with increasing reflects the fact that, for (3 = 0, LLM has a 
significantly stronger effect on Asp than on the Asr, and so tends to stabilize the spin- 
polarized excitation. 

For the more experimentally relevant case (3 = 1.5 Iq the effect of LLM on Asp and Asr 
is much weaker than for j3 = and the difference in the two gap energies, in units of e^/elo, 
is roughly constant. Again using GaAs parameters we find i?c — 4 T for = 2 and -Be — 3 
T for Ts = 10. This result for the crossover field is in reasonable agreement with previous 
calculations which included the thickness correction but not LLM . The new result here 
is that, when thickness is included. Be is only weakly dependent on LLM. 

Figure 3 shows mixed estimates of the density profiles of the spin-polarized and spin- 
reversed excited states at z/ = 1/3 as a function of 6 for = 1 and = 20. The results are 
for = 20 and /5 = 0. Note that for = 20 the quasielectron and quasihole induce ripples 
in the density. This is a consequence of the increased Wigner crystal-like correlations in the 
FQHE fluid induced by LLM. 

The effect of LLM on Asp and Asr can be understood qualitatively using Fig. 3 as 
follows. For the spin-polarized state the quasielectron charge is concentrated in a ring 
centered at the top of the sphere. As increases this ring of charge spreads out and the 
dip at 6* = deepens. This evolution of the quasielectron charge occurs because, as higher 
Landau levels are mixed into the wave function, the charge is free to spread out, lowering 
its Coulomb energy at the cost of some kinetic energy. For the spin-reversed quasielectron 
the charge is spread out more uniformly and the Coulomb energy is less than that of the 
spin-polarized quasielectron. There is therefore less of a tendency for the spin-reversed 
quasielectron charge to spread out with increasing r^, and so this excitation is less affected 
by LLM. The weakening of the effect of LLM on the energy gaps due to the thickness 
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correction can also be understood along similar lines. For [3 = 1.5 Iq the short-range part 
of the Coulomb interaction is softened, and the potential energy of both the spin-reversed 
and spin-polarized quasielectrons are reduced. There is therefore less of a tendency for the 
charge of these excitations to spread out with increasing r^, and so, again, the effect of LLM 
is weakened. 

For typical carrier densities the magnetic field at z/ = 1/3 is greater than Be and the 



transport gap is set by Asp. This gap has been measured in both n-type ||T6| and p-type ||T 
GaAs quantum wells, with typical results, for the highest quality samples, of An — 0.05e^/e/o 
and Ap ~ 0.023e^/e/o, respectively. The factor of 2 reduction of the energy gap from n- 



type {vg ~ 2) to p-type {vg ~ 10) samples has been attributed to the increased LLM [117 



However, our results show that LLM has only a weak effect on Asp when the thickness effect 



is included, in agreement with previous calculations [0,|1^,0. Thus, while our energy gap 
for ^ 2 is close to the experimental value, our result for ^ 10 is off by roughly a factor 
of 2. This discrepancy between theory and experiment is most likely due to disorder, the 
effects of which on the energy gap are still poorly understood. 

To summarize, a generalized DMC method for solving the many-body Schrodinger equa- 
tion on curved manifolds has been introduced and used to perform a 'fixed-phase' simulation 
of the FQHE on the Haldane sphere. The effect of LLM on the z/ = 1/3 energy gap, and the 
relative stability of the spin-polarized and spin-reversed quasielectron states have been inves- 
tigated using the new method. We believe that the generalization of the 'fixed-phase' DMC 
method to the spherical geometry presented here will be useful for many future numerical 
studies of the FQHE. 
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FIG. 1. Ground state energy per particle for = 1/3 as a function of r^. The dashed hne 



is the variational result of Price et al. [13| and the solid line is a least square fit of second degree 
polynomial in to our fixed-phase DMG results for = 1, 5, 10, 15 and 20 (dots). Statistical 
errors are smaller than symbol sizes. 
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FIG. 2. v = 1/3 energy gaps for creating a quasielectron — quasihole pair at opposite poles of 
the sphere vs. Vg- Results are given for both a spin-polarized (dots) and spin-reversed (diamonds) 
quasielectron for thickness parameter (a) /3 = 0, and (b) (3 = 1.5 Iq. Results are for 20 electrons. 
The lines are guides to the eye. 
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FIG. 3. Mixed estimates of the density profiles of u = 1/3 excited state wave functions on the 
sphere with a quasielectron at the top of the sphere {9 = 0) and a quasihole on the bottom of the 
sphere {6 = tt) with = 1 (dashed hne) and = 20 (sohd hne), for (a) the spin-polarized excited 
state, and (b) the spin-reversed excited state. Results are for 20 electrons. 
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